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Abstract 

We present preliminary results of the analysis of B — ► D°h decays, with h = ir,K and the D° 
reconstructed in the channels K~ir + , K~Tr + ir + ir~ , K~ir + ir° and in the CP eigenstate K~K + , using 
data collected by the BABAR detector during the years 2000-2002 at the PEP-II asymmetric-energy 
B Factory at SLAC. We have measured the ratio of the branching fractions 
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and the direct CP asymmetry 
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1 Introduction 



During the last ten years there has been growing theoretical interest in knowing the decay rates 
for the processes B~ -» D°K~ Q and B ± -> D^pK*, where D° CP indicates the CP -even or CP- 
odd states (D° ± D°)/\^2. These modes are key ingredients for some of the recently proposed 
methods for extracting the angle 7 of the Cabibbo-Kobayashi-Maskawa quark mixing matrix in a 
theoretically clean way gj. 

In this paper we present an analysis of the decay B~ — > D°K~ in which the D° is reconstructed 
in the non-CP eigenstates K~tt + , K~tt + tt + -k~ , K~7t + -k°, or in the CP-even eigenstate K~ K + . The 
ratio R between the branching fractions of B~ — > D°K~ and B~ — ► D°ir~ 

_ B(B~ -» Z?Qg-) 

S(B~ -► £>%-) 1 j 

is measured, along with the ratio Bcp between the branching fractions of B — > -Dgp if and 

_ fi(B~ -> ^ P K-) + B(B+ -> 
CP _ B(B~ ^ P vr-) + £(B+ -» ^p7r+) ' 1 j 

The yields of B~ — > Dq P K~ and B + — > -D^p -ftT + are separately extracted, and the CP asymmetry 

= g(B~ -» D°p K-j - B(B+ -» D°p if +) 
CP e(B-^ J D°p^-)+e(P+^ J D°pK+) ' 



is measured. 

The B~ — > D°K~ decay was first observed by the CLEO collaboration [||. Using a sample 
of 3.1 fb" 1 collected at the T{4S) resonance, CLEO measured R = (5.5 ± 1.4 ± 0.5)%. Recently, 
based on 10.4 ftT 1 collected at the T(4S), Belle has measured R = (7.9 ±0.9 ±0.6)% (§. Belle also 
reported, using a sample of 29.1 fb _1 , measurements of the direct CP asymmetries for the CP-even 
and CP-odd modes, A C p(+l) = (0.29 ± 0.26 ± 0.05)% and A C p{-l) = (-0.22 ± 0.24 ± 0.04)% g. 



2 The BABAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring during the years 2000-2002. The sample corresponds to an integrated luminosity of about 
75 fb -1 accumulated at the T(4S) resonance ( "on-resonance" ) and about 10 fb -1 accumulated 
at a center-of-mass (CM) energy about 40 MeV below the T(4S) resonance ("off-resonance"). 
Data taken below the Y(4S) are used for continuum background studies. The on-resonance sample 
corresponds to (81.1 ±0.9) x 10 6 BB pairs. The results regarding the non-CP modes are based on a 
subsample of data collected in the years 2000-2001; this corresponds to an integrated luminosity of 
about 56 fb -1 accumulated at the T(4S) resonance, which in turn corresponds to (61.2 ±0.7) x 10 6 
BB pairs. 

PEP-II is an e + e~ storage ring operated with asymmetric beam energies, producing a boosted 
(/?7 = 0.55) T(4S) along the collision axis. BABAR is a solenoidal detector optimized for the 
asymmetric beam configuration at PEP-II and is described in detail in Ref. ||. Charged particle 
(track) momenta are measured in a tracking system consisting of a 5-layer, double-sided, silicon 

1 Charge conjugation is implied here and throughout this paper unless explicitly stated. 



8 



vertex tracker (SVT) and a 40-layer drift chamber (DCH) filled with a gas mixture of helium and 
isobutane, both operating within a 1.5 T superconducting solenoidal magnet. Photon candidates are 
selected as local maxima of deposited energy in an electromagnetic calorimeter (EMC) consisting of 
6580 CsI(Tl) crystals arranged in barrel and forward endcap subdetectors. Particle identification is 
performed by combining information from ionization measurements (dE/dx) in the SVT and DCH, 
and the Cherenkov angle 9c measured by a detector of internally reflected Cherenkov light (DIRC). 
The DIRC system is a unique type of Cherenkov detector that relies on total internal reflection 
within the radiating volumes (quartz bars) to deliver the Cherenkov light outside the tracking and 
magnetic volumes, where the Cherenkov ring is imaged by an array of ~ 11000 photomultiplier 
tubes. 

3 Event selection 

We reconstruct B mesons decaying to a D° meson and a charged prompt track h, where h is 
a pion or a kaon. D° meson candidates are reconstructed in four decay modes: D° — » K~7r + , 
D° -> K-TT+TT+7T-, D° -> K-it + it and D° -> K~K+. The decay B~ -» D°tt- is used to evaluate 
detector resolutions, systematic uncertainties and to normalize branching fractions. 

All charged tracks are reconstructed in the drift chamber and/or the vertex detector, and 
their parameters determined with the pion mass hypothesis. In order to reduce the combinato- 
rial background, only charged tracks with momentum greater than 150 MeV/c are used in the 
reconstruction of D° — > K~TT + ir + ir~ and D° — > K~tt + tt°; the prompt track h is required to have 
momentum greater than 1.4 GeV/c. Particle ID information from the drift chamber (dE/dx) and, 
when available, from the DIRC is required to be consistent with the kaon hypothesis for the K 
meson candidate from the D . In order for the prompt track h to be identified as a pion or a kaon 
we require that its Cherenkov angle be reconstructed in the DIRC with at least five photons. We 
reject candidate tracks whose Cherenkov angle is within 3 standard deviations (a, where a means 
the experimental resolution) from the expected angle for the proton hypothesis, and candidate 
tracks that are identified as electrons by the DCH and the EMC. 

Candidate 7r° mesons are reconstructed from a combination of two photon candidates. Photon 
candidates are selected as showers in the EMC that have the expected lateral shape, are not matched 
with any charged track, and have a minimum energy of 70 MeV. The 77 invariant mass is required 
to be in the range 124-144 MeV/c 2 , and the total energy of the 7 pair must exceed 200 MeV. 
The 7T° — * 77 candidates are then kinematically fit with their mass constrained to the nominal tt° 
mass J?]]. 

The invariant mass of D° — ► K~ir + , D° — ► K~ir + ir + 7r~ and D° — > K~ K + candidates is 
required to be within 3 a of the mean mass value. As the combinatorial background of the D° — > 
K~it + it decays is larger due to the presence of a neutral pion, a 2-a cut on the invariant mass is 
applied for this mode. D° candidates are then kinematically fit with their mass constrained to the 
nominal D° mass Q, in order to improve momentum determination. 

We reconstruct B meson candidates by combining a D° candidate with a track h. For the 
non-CP modes, the charge of the bachelor track h must match that of the kaon from the D° meson 
decay. We exploit our knowledge of the initial state momentum and energy to select B meson 
candidates by defining two largely uncorrelated variables. The first is the beam- energy substituted 
mass tties = \J ~ Pb > where E^ = (s/2 + p, • p# )/!?,, yfs and E$ are the total energies of the 
e + e~ system in the CM and lab frames, respectively, and p, and p# are the momentum vectors 
in the lab frame of the e + e~ system and the B candidate, respectively. The second is the energy 



9 



difference AE, which is defined as the difference between the energy of the B candidate and half 
the energy of the e + e~ system, computed in the CM system. The toes resolution is dominated 
by the beam energy spread, while for AE the main contribution comes from the measurement of 
particle energies in the detector. 

The toes distribution for B~ —* D°h~ signals does not depend on whether the prompt track 
h is a pion or a kaon, nor, to first order, on the D° momentum resolution. This has been checked 
both for signal Monte Carlo and for data. The m^s resolution for the B~~ — > D°h~ channels is 
found to be 2.6MeV/c 2 , regardless of the D° decay mode. 

In contrast, the AE distributions depend on the nature of the bachelor track h and on the 
resolution for the D° meson momentum. We evaluate AE with the kaon mass hypothesis (and 
denote it by AEk to avoid confusion) in such a way that the distributions are centered near for 
B~ — > D°K~ events and shifted by approximately 42MeV for B~ — ► D -k~ events. The AEk 
resolution is typically 17MeV for all D° decay modes. 

B candidates are selected in the range 5.2 < m ES < 5.3GeV/c 2 and -0.100 < AE < 0.130 GeV. 
For events with multiple B candidates, the best candidate is chosen based on the values of the D° 
invariant mass and toes- 

4 Background rejection 

The physics backgrounds for the considered modes originate both from the continuum production 
of light quarks, e + e~ — > qq (q = u,d,s,c), and from BB events. In the center-of-mass frame, the 
continuum background typically exhibits a two-jet structure; in contrast, the low momentum of the 
B meson in the decay T(4S) — > BB leads to a more spherically symmetric event. We exploit this 
topology difference between signal and continuum background by making use of two event-shape 
quantities. 

The first variable is the normalized second Fox- Wolfram moment ||, R2 = H*, where Hi is the 
£-order Fox- Wolfram moment. R2 is required to be less than 0.5 for all the selected events. 

The second quantity is the angle 6t between the thrust axes, evaluated in the center-of-mass 
frame, of the B candidate and the remaining charged and neutral particles in the event. The 
absolute value of the cosine of this angle is strongly peaked near 1 for continuum events and is 
approximately uniform for BB events. | cos 6t\ is required to be less than 0.9 for the D° — > K~ir + 
mode, and this value is tightened to 0.7 for D° — > K~it + it + it~ and D° — > K~it + it modes, which 
suffer from larger combinatorial background. 

For the D° —* K~K + mode an additional quantity, the D° rest frame decay angle 9kk, is used 
in conjunction with | cos 6t\- The angle 0k k is defined as the angle between the direction of the 
D° calculated in the rest frame of the B and the direction of one of the decay products of the D° 
calculated in the rest frame of the D°. The distribution of cosOkk is flat for signal and peaked at 
±1 for fake D° background. | cos 9t\ and cosOkk are uncorrelated for signal but not for continuum 
background. This correlation is exploited to make a more efficient cut in the cos 9t — cos Ork plane. 

The main contributions from BB background come from the processes B~ — > D*°h~ (h = 
7r,K), B~ — > D°p~ and from mis-reconstructed B~ — > D°h~ decays. Background from the 
charmless three-body process B~ — > K~ K + K~ is potentially the most critical for the B~ — > 
D°K~, D° — > K~ K + decay, because it consists of three kaons coming from a B meson and 
hence is characterized by the same AE and ttt-es distribution as the signal. Results from recent 
studies of the resonant composition of such a decay have been used to estimate this non-negligible 
background. 
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The total number of B candidates and the final selection efficiency e for each mode are sum- 
marized in Table p]. The small differences between B~ — > D Q ir~ and B~ —* D°K~ efHciences 
arise both from reconstruction (kaons have lower dE/dx than pions and hence fewer hits associated 
to the track) and from candidate selection (the rejection of electrons and protons has different 
efficiencies for kaons and pions). 

Table 1: Number of candidates selected for the maximum- likelihood fit and final selection efficien- 
cies. The selection efficiency has been evaluated on simulated events. 





decay mode 


events selected 


e{B~ -» D u 7r~) 


e(B~ -» D^K-) 


fCdt (fb^ 1 ) 


K- 


7T+ 


12606 


(43.65 ±0.12)% 


(42.17 ±0.31)% 


56.0 


K- 




9782 


(14.53 ±0.17)% 


(13.55 ± 0.22)% 


56.0 


K- 


7T+7T 


8177 


(8.04 ±0.11)% 


(8.12 ±0.18)% 


56.0 


K- 


K+ 


2389 


(34.5 ± 0.2)% 


(33.7 ±0.3)% 


74.8 



5 Signal extraction 

For each D° decay mode an extended unbinned maximum-likelihood fit to the selected data events 
determines the signal and background yields rtj (i = 1 to M, where M is the total number of signal 
and background species). We consider two kinds of signal events, B~ — > D°tt~ and B~ — ► D°K~ , 
and four kinds of background events: candidates selected either from continuum or from BB events, 
in which the prompt track h is either a pion or a kaon. In the case of B~ — > D°h~ , D° — > K~K + 
events, the fit is also performed on the B + and B~ subsamples separately. 

The input variables to the fit are toes) AEr, and a particle identification probability for the 
bachelor track h based on the Cherenkov angle 6c, the momentum p at the DIRC and the polar 
angle of the track. The extended likelihood function C is defined as 



C = exp 




■ M 
i=l 



(4) 



The exponential factor in the likelihood accounts for Poisson fluctuations in the total number of 
observed events N. The M functions Vi(xj; &i) are the probability density functions (PDFs) for the 
variables Xj, given the set of parameters 5j. They are evaluated as a product V% = Vi(AEk, mEs) x 
"Pi (DIRC), since m^s and AEx are not correlated with the Cherenkov angle of the prompt track. 

The correlation between toes an( l AEk can be neglected for signal events and for the com- 
binatorial background from the continuum: V^AE^, Wes) = 'Pi(AEx) x V^m^s). However it 
cannot be neglected for the BB background component, for which we use a two-dimensional PDF 
determined from simulated events through a method based on the Kernel Estimation [|l0| technique. 

The parameters for the AEx and tt^es distributions for continuum background are determined 
from off-resonance data. The background shape in AEk is parameterized as a linear polyno- 
mial, while the ?tt,es shape is parameterized by an ARGUS threshold function |[l]] /(mEs) (X 
Jtt-esVI — x 2 exp[— £(1 — x 2 )], where x = mEs/^o a n <4 tuq is the average CM beam energy. 

The tties PDF for signal events is determined from a pure sample of B~ — > D°ir~ , D° — ► K~ir + 
decays selected from on-resonance data by applying the standard selection plus a 3 o cut on AE W . 
BB background is negligible in this channel; therefore the tties distribution is fit by a Gaussian 



11 



signal over a background described by an ARGUS function. The resulting Gaussian shape is used 
to parameterize the tties PDF for the B~ — ► D°h~ signals. 

The parameterization of the AEk distribution for B — ► D°K~events is deduced from that 
of AE n for B — > D°7r~ events in the data. The AE^- distribution for a pure sample of B~ — > 
.D 7r _ decays, selected by applying the standard selection plus a 3 cr cut on tbes and the requirement 
that the bachelor track h not be consistent with the kaon hypothesis, is fit by a Gaussian signal 
over a linear background. The resulting Gaussian shape is used to parameterize the AEk PDF 
for the B~ — > D°K~ signal. We also use the AE W resolution together with the known momentum 
spectrum of the bachelor track in the laboratory frame to generate a "translated" AEk distribution 
for B — > D°tt~ events: for each generated value of AE W and p, we calculate the energy shift 
according to Ai^hift = AEk — AE n = 7 {^\Jrn K + p 2 — \/m 2 + p 2 ^j , where 7 is the Lorentz boost 
of the center-of-mass frame. The resulting AEk distribution is empirically parameterized by the 
sum of two Gaussians. 

The 6c PDFs are derived from kaon and pion tracks in the momentum range of interest from 
a sample of D* + D°tt + (D° R-tt + ) decays. Since the 0c resolution depends in principle 
uniquely on the number of found photons, and the number of expected photons depends on both 
the track polar angle and momentum, we use this sample to parameterize the 6c resolution as a 
function of track polar angle and momentum extrapolated to the DIRC. 

The results of the fit are summarized in Table || and the AEk distributions are shown in Fig. |l| 
for all the D° modes. The projection of the likelihood fit is overlaid on the distributions. The 
individual contributions of each of the signals and the total background are shown. 

In order to increase the relative fraction of signal B — > D°K~ events for illustration only, 
candidates are selected with the further requirements that the prompt track be identified as a kaon 
and |mES — ( m Es)| < 3<r. The AEk distributions of the selected candidates are shown in Figure |2| 
for all modes. The peak around zero from B — > D°K~ candidates is visible, but the efficiency for 
the signal is decreased by 13% with respect to the standard selection outlined in Sees. ^ and ^. 

Table 2: Results from the maximum-likelihood fit. For the D° — > i^-fT+mode we quote the results 
for the fits performed on the whole sample and on the B + and B~ subsamples. 



D [) decay mode 


N(B~ -»■ D u tt-) 


N(B~ D°K-) 


JCdt (fb- 1 ) 


K-TT+ 


4440 ± 69 


360 ± 21 


56.0 




2914 ± 56 


242 ± 18 


56.0 


K-1T+1T 


2650 ± 56 


208 ± 18 


56.0 


K-K+ 


508 ± 24 


37 ±8 


74.8 


K-K+ [B+] 


254 ± 17 


15 ±6 


74.8 


K-K+ [B~] 


254 ± 17 


22 ±6 


74.8 



6 Systematic studies 

Systematic uncertainties on the ratio R and on the CP asymmetry Acp arise primarily from un- 
certainties on rii due to imperfect knowledge of the PDF shapes. In this section we describe how 
the different contributions to these uncertainties have been evaluated; the results are summarized 
in Tables | and |. 
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Figure 1: AEk distribution of the B — ► D°h candidates selected in the data sample. As 
described in the text, tracks h consistent with the electron or proton hypothesis have been explicitly 
removed. Top left: D° -> K~7r + ; top right: D° -> K--k + it + it-; bottom left: D° -> if-vr+vr ; 
bottom right: D° — ► K~ K + . The solid curve represents the projection on the AEk axis of the 
resulting fit probability density function, scaled by the number of candidates in the sample. The 
contribution from the two signals and from the total background is also shown. 



In the case of the mss an d AEk PDFs for the different signal and continuum background 
contributions, we vary the PDF parameters — which have been determined through a fit to AEk 
and mES distributions in data or Monte Carlo events — by one statistical error. 

The systematic uncertainty from the BB background must be treated as a special case, because 
the parameterization is performed through the Kernel Estimation and not through an analytical 
function. To study the systematic uncertainty associated with the uncertainty on the shape of the 
BB background we generate 500 different samples of BB events and for each one the correspond- 
ing PDF is used in the maximum-likelihood fit. The width of the distribution of the difference 
between the new yields and the original yields is used as an estimate of the systematic uncertainty. 
Uncertainties arising from the imperfect knowledge on the branching fractions of the different chan- 
nels that contribute to the BB background are also taken into account (this is important for the 
D° — ► K~K + mode, where one of the main sources of uncertainty is the expected number of 
B — > KKK background events). 

The parameterization of the particle identification PDF is performed by fitting with a Gaussian 
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AE K (GeV) AE K (GeV) 

Figure 2: AEk distribution of the B~ — > D°h~ candidates selected in the data sample by requiring 
that the prompt track h be identified as a kaon and |mES — ("1es)| < 3<r. Top left: D° — ► K~ir + ; 
top right: L>° -> K-vr+vr+vr-; bottom left: D° -> if-vr+Tr ; bottom right: L>° -> R-R+. The 
peak around zero from 5~ — > D°K~ candidates is visible. 



function the background-subtracted distribution of the difference between the reconstructed and 
expected Cherenkov angle of the charged tracks from D° decays in the D* + — > D°tt + (D° — > K~ir + ) 
control sample, in bins of momentum and polar angle. Therefore, the significant parameters for 
each momentum-polar angle bin are the mean and the width of the fitted Gaussian. We estimate 
the systematic error associated with the particle identification PDF by varying by one statistical 
error the Gaussian parameters of each bin, while all the others are kept fixed at their central value. 

The various sources of systematic errors are assumed to be uncorrelated. The total systematic 
error is obtained by summing in quadrature the individual contributions. 

Systematic uncertainties on the charge asymmetries are added in quadrature with the limit 
on intrinsic charge bias in the detector (0.05), measured on data for the processes B~ — > D°n~ 
[D° K~tt + ], B- -» D°K~ [D° -► R-7T+] and B~ D°ir~ [D° K~K+], where the CP 
asymmetry is expected to be zero. 

7 Results 
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Table 3: Systematic uncertainties on the ratio R of the branching fractions of B — > D°K and 
B-^D°7T-. 



Parameter 


Uncertainty on R (%) 


AE K (B- -> D°h~) 


±0.13 


m ES (B- -> D°h-) 


±0.01 


qq background AEk 


±0.01 


qq background rngs 


±0.01 


BB background AEk vs mES 


±0.04 


Particle identification 


±0.16 


Background composition 


±0.03 


Total 


±0.20 



Table 4: Systematic uncertainties on the CP asymmetry in the decay B — ► D K , with D' 
K~K+. 



Parameter 


Uncertainty on Aqp (%) 


AE K (B~ -> D°h~) 


+0.5 
-0.0 


m ES (B- -> D°h-) 


+0.5 
-0.1 


qq background AEk 


+0.3 
-0.0 


qq background th-es 


+0.1 
-0.7 


BB background AEk vs tues 


±0.3 


Particle identification 


+4.1 
-0.9 


Background composition 


+6.2 
-4.4 


Detector charge asymmetry 


±5.0 


Total 


+9.0 
-6.8 



7.1 Measurement of the ratio B(B~ -> D°K-)/B(B- -> D°n-) 

The ratio of the branching fractions of the decays B~ — > D°n~ and — > D°K~~ is calculated 
separately for the three non-CP Z)° decay channels and for the D° — > K~ K + mode. The ratio 
is computed from the number of B~ — > D°K~ and B _ — > D°ir~ mesons estimated with the 
maximum-likelihood fit and listed in Table ^. The resulting ratios are scaled by a correction factor 
that account for small differences in the efficiency between B~ — > D°K~ and B~ — > D°tt~ selection, 
estimated with signal Monte Carlo samples. The results for the ratio of the branching fraction of 
the decay B~ — ► D°tt~ and — > D°K~ with Z)° decaying to the non-CP modes, measured 
with 56.0 fb _ , are listed in Table ||; they represent three independent measurements of the same 
quantity, and they are statistically consistent. The weighted mean of the three measurements gives 
the result for the non-CP modes: 

R = (8.31 ± 0.35 ± 0.20)% . (5) 
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The resulting ratio of the branching fractions of the D° — > K K + mode, measured with 74.8 fb , 
is also reported in Table || and is found to be: 

R CP = (7.4 ± 1.7 ± 0.6)% . (6) 

In the evaluation of the statistical error on R the correlation between the numbers of B~ — ► D°K~ 
and B — > D ir~ events obtained from the fit has been taken into account. 

Table 5: Measured ratio B(B~ — > D°K~)/B(B~ — ► D°7r~) for different Z?° decay modes. 



— > Z? u /i decay mode 


ratio 


fCdt (fb^ 1 ) 


-» K~tt + 


(8.4 ±0.5 ±0.2) % 


56.0 


D° -> i^-vr+vr+vr- 


(8.7 ±0.7 ±0.2)% 


56.0 




(7.7 ±0.7 ±0.3)% 


56.0 


weighted mean 


(8.31 ±0.35 ±0.20) % 


56.0 




(7.4 ±1.7 ±0.6)% 


74.8 



7.2 Measurement of the direct CP asymmetry 

The direct CP asymmetry for the B~ —* D°K~, D — > K~K + decay is calculated from the 
measured yields of positive and negative decays reported in Table |2[ The resulting asymmetry, 
measured with 74.8 fb -1 , is 

A CP = 0.17 ± 0.23 ±°$ . (7) 

We checked for detector charge bias by measuring the asymmetries for the processes B — ► D°ir~ 
[D° -» K-7T+], 5- D°K~ [D° -» and S" -> P> vr- [D° -> where the 

CP asymmetry is expected to be zero. The measured values, (—2.4 ± 1.3)%, (—0.6 ± 5.0)% and 
(0.0 ± 4.7)%, are consistent with zero. 



8 Summary 

The B~ — ► D°K~ decays with £) selected in the channels D° — > K~tt + , D° — > Er _ 7r + 7r + 7r _ , 
-7T + 7T have been reconstructed on a data sample of 56.0 fb . The ratio R of the 
branching fractions B(B~ — » D°K~) and B(B~ —* D°7r~) has been measured to be 

0fB~ _, D°K-) 
R = A_ _ = ( 8 .3l ± 0.35 ± 0.20)% . 

The -B~ — > D°K~ decays with D° decaying to the CP-even eigenstate D° — ► K~K + have been 
reconstructed on a data sample of 74.8 fb -1 . The yield has been measured separately for positive 
and negative B mesons. The total yield is 

N± = N B± ^ D o cpK± = 36.8 ± 8.4 ± 4.0 . 
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The ratio Rcp of the branching fractions B{B^ — ► D^pK^) and B(B^ — ► D^pir^) has been 
measured: 



7? - B(B~ -» D° P ^-) + (ff+ -» £° P A-+) _ 
^ CP = B(S-^ J D^-) + B(S+-2?^ " (7 - 4 ± L7 ± °- 6) ' 
The direct CP asymmetry has been measured: 



o 



AcP = B(B~ - D° CP K~) + B(B+ - D^K+) ~ ^ ± ^ -°' 07 ' 
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